Emergence of large scale structure 
in planetary turbulence 
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Turbulent flows are observed to self-organize into large scale structures such as zonal jets and 
coherent vortices. The simplest model that retains the relevant dynamics is a barotropic flow in 
a /3-plane channel with turbulence sustained by random stirring. Non-linear integrations of this 
model show that as the energy input rate of the forcing is increased, the homogeneity of the flow is 
broken first by the emergence of non-zonal, coherent, westward propagating structures and second 
for larger energy input rates by the emergence of zonal jets. We use a non-equilibrium statistical 
theory, Stochastic Structural Stability Theory (SSST), to address the emergence of these coherent 
structures. Employing the tools of SSST, we construct a coupled dynamical system governing the 
joint evolution of the coherent flow and of the second order statistics of the turbulent eddies. We 
then treat the structural stability of a homogeneous equilibrium with no mean flow analytically. 
We find that non-zonal and zonal coherent structures appear as the result of structural instability 
and equilibrate at finite amplitude. We also find that SSST accurately predicts the critical energy 
input rates for the emergence of both non-zonal and zonal coherent structures in the non-linear 
integrations as well as the characteristics (scale, amplitude and phase speed) of these structures. 
PACS numbers: 47.27.eb, 47.20.Ky, 47.27.De, 52.35.Mw, 89.75.Fb, 92.60. Bh, 92.10.A- 



Turbulcnce in planetary atmospheres and in plasma 
flows is commonly observed to be organized into large 
scale unidirectional (zonal) jets with long-lasting coher- 
ent eddies or vortices embedded in them [1, 2]. The jets 
control the transports of heat and chemical species in 
planetary atmospheres and separate the high tempera- 
ture plasma from the cold containment vessel wall in 
magnetic plasma confinement devices. It is therefore 
important to understand the mechanisms for the emer- 
gence, equilibration and maintenance of these coherent 
structures. In this Letter we present a theory for the 
mean turbulent state that is able to predict the regime 
changes occurring in the turbulent flow as well as the am- 
plitude, structure and propagation characteristics of both 
the zonal jets and the non-zonal coherent structures that 
form in the flow. We then test this theory against non- 
linear simulations in a simple model of forced planetary 
and plasma turbulence. 

The simplest model that captures the turbulent dy- 
namics and its interaction with the zonal jets and the co- 
herent structures, is the stochastically forced barotropic 
vorticity equation on a plane tangent to the surface of a 
rotating planet: 

dtC + - ^yC, + * - -r( - f A 2 C + /■ (1) 



meridional (north-south) direction, j3 — 2ilcos(j>o/R is 
the gradient of planetary vorticity, Q the rotation rate 
of the planet, </>q the latitude of the /3-plane and R is 
the radius of the planet. Equation (1) governs the dy- 
namics of non-divergent motions at the midlatitudes of 
the planet. It is also the infinite effective Larmor ra- 
dius limit of the Charney-Hasegawa-Mima equation that 
governs drift- wave turbulence in plasma flows. We are 
assuming linear damping with coefficient, r, representing 
the Ekman drag induced by the horizontal boundaries 
and hyper-diffusion with coefficient v that dissipates the 
energy flowing into unresolved scales. The forcing term / 
is necessary to sustain turbulence, and may parameterize 
processes that have not been included in the dynamics, 
such as the vorticity forcing from small scale convection. 
In many previous studies, this exogenous excitation is 
taken as a temporally delta correlated and spatially ho- 
mogeneous and isotropic random stirring. We will follow 
the same forcing protocol in this Letter and consider an 
isotropic ring forcing that is injecting energy at rate e in 
a narrow ring of wavenumbers of width AKf around the 
total wavenumber Kf. 

Equation (1) is solved in a doubly periodic domain 
of size 2ir x 2ir. The calculations presented in this Let- 
ter are for /3 = 10, r = 0.01, v = 2 • 10~ 6 , K f 



The relative vorticity is ( = Aip, where vb is the stream- and AKf = L Thc resu i ts discussed were verified 



function, A 



dyy is the horizontal Laplacian, 
x is in the zonal (east-west) direction and y is in the 
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be insensitive to the parameters and the forcing proto- 
col chosen. To illustrate some of the characteristics of 
the turbulent flow and the emergence of coherent struc- 
tures, we consider two indices. The first is the zonal mean 
flow index [3] defined as thc ratio of the energy of zonal 



jets over the total energy, zmf = 
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FIG. 1. The zmf (red lines) and nzmf (blue lines) indices de- 
fined in the text as a function of energy input rate e/e c for the 
non-linear (solid lines) and SSST (dashed lines) integrations. 
The critical value e c is the energy input rate at which the 
SSST predicts structural instability of the homogeneous tur- 
bulent state. Zonal jets emerge for e > e n i, with e n i = 15.9e c 
for the parameters in the simulation. 
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FIG. 2. Time averaged energy power spectra, log 10 (-E(fc, /)), 
obtained from non-linear simulation of Eq. (1) at (a) e/e c — 
2.6 and (b) e/e c = 30. In (a) the flow is dominated by a 
(\k\, \l\) — (1, 5) non-zonal coherent structure. In (b) the flow 
is dominated by a coherent zonal flow at (k, \l\) = (0, 3). 



is the time averaged energy power spectrum of the flow 
and k, I are the zonal and meridional wavenumbers re- 
spectively. The second is the non-zonal mean flow index 
defined as the ratio of the energy of the non-zonal modes 
with scales lower than the scale of the forcing over the 

-E;S(fe=o,o) 

-• Fig- 
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total energy: nzml = — w 

ure 1 shows both indices as a function of the energy input 
rate e. 

For e smaller than a critical value e c , the turbulent 
flow is homogeneous and remains translationally invari- 
ant in both directions. When e > e c , the translational 
symmetry of the flow is broken and non-zonal structures 
form with scales larger than the scale of the forcing. The 
time averaged power spectrum, shown in Fig. 2(a) for 
e = 2.6e c , has pronounced peaks at (k, I) = (±1, ±5) that 
correspond to coherent structures propagating westward 
(cf. Fig. 3(a),(b)) with approximately the Rossby wave 
phase speed for this wave. However, at larger e the prop- 
agation speed of these structures departs from that of 
Rossby waves. The presence and properties of such non- 




FIG. 3. (a) Snapshot of the streamfunction tp(x,y,t) and 
(b) Hovmoller diagram of ip(x,y = ir/A,t) obtained from non- 
linear simulation for e/e c = 2.6. The thick lines in (b) cor- 
respond to the phase speed obtained from the stability equa- 
tion (6).(c) Snapshot of the streamfunction ijj(x,y,t) and (d) 
Hovmoller diagram of the z-averaged ip(y,t), which is ob- 
tained from non-linear simulations at e/e c = 30. 



linear waves in similar simulations, were also reported 
recently [4] . For e > e n i indicated in Fig. 1 , robust zonal 
jets emerge. For example the peaks at (k, I) = (0, ±3) in 
the spectrum of Fig. 2(b) correspond to coherent zonal 
jets (cf. Fig. 3(c), (d)). From Fig. 1 we see that while 
the jets contain over half of the total energy, substantial 
power remains in non-zonal structures. Previous studies, 
refer to the coherent non-zonal structures in this regime 
in which strong zonal jets dominate the flow (referred 
to as zonostrophic regime [5]) as satellite modes [6] or 
zonons [4]. 

The emergence of the jets has been described in pre- 
vious studies in terms of an anisotropic inverse energy 
cascade [7-9], or in terms of inhomogeneous mixing of 
vorticity [10], or in terms of a direct transfer of energy 
from small scale waves into the zonal jets, through cither 
non-linear interactions between finite amplitude Rossby 
waves [11, 12], or through shear straining of the small 
scale waves by the jet [13]. In addition, the mechanism 
for the emergence of the non-zonal structures remains elu- 
sive. Statistical equilibrium theory applied in the absence 
of forcing and dissipation, has been able to predict both 
jets and coherent vortices as maximum entropy struc- 
tures [14] and a recent study has shown correspondence of 
the theoretical results with non-linear simulations in the 
limit of weak forcing and dissipation [15]. However, the 
relevance of these results in planetary and plasma flows 
that are strongly forced and dissipated and are therefore 
out of equilibrium remains to be shown. In this Letter 
we present results of an alternative, non-equilibrium sta- 
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tistical theory, that is termed as Stochastic Structural 
Stability Theory (SSST) [16, 17] or Second Order Cumu- 
lant Expansion theory (CE2) [3, 18, 19]. While recent 
studies have demonstrated that SSST can predict the 
structure of zonal flows in turbulent fluids [3, 20, 21], 
the results presented in this Letter, demonstrate that an 
extended version of SSST can predict the emergence of 
both zonal and non-zonal coherent structures in plane- 
tary turbulence and can capture their finite amplitude 
manifestations. 

SSST describes the statistical dynamics of the first two 
equal time cumulants of Eq. (1). The first cumulant is 
Z(x, t) = (C) (the brackets denote an ensemble average) 
and the second cumulant C(xi,x 2 ,i) = (C1C2) i s a func- 
tion of the vorticity deviation Q = (i — Zi at the two 
points Xi = (x i: yi) (i = 1,2) at time t. It can be read- 
ily shown that the equations for the evolution of the two 
cumulants are: 

d t Z + UZ X + V{fi + Z y ) +rZ + vA 2 Z = 

= d x (d yi A^C) xi=x2 - d y (d Xl A^C) xi=x2 , (2) 
d t C= (A 1 +A 2 )C + E. (3) 

The linear operator 

Ai = -Uid Xi - V t d Vi -(J3 + Z yi )d Xi A7 1 + 

+ Z Xi d yi A7 1 -r-vAl (4) 

acts at the points Xj = (xi,yi) and governs the dynam- 
ics of linear perturbations about the instantaneous mean 
flow U = [U,V] = [-dy (tp),d x (tp)]. In (3), S contains 
the covariance of the external forcing and terms related 
to third order cumulants. A second order closure is ob- 
tained if the third order cumulant is ignored and 5 is set 
to be the spatial covariance of the stochastic forcing /. 
In most earlier studies of SSST or CE2, the ensemble av- 
erage was assumed to represent a zonal average. In this 
Letter, we adopt the more general interpretation that the 
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where K 2 = k 2 + l 2 , K 2 = (k + n) 2 + (I + to) 2 , N 2 = n 2 + 
to 2 , k + = k + n/2 and l + = I + to/2 and the summation 
is over integers (k, I) satisfying \K — Kf \ < AKf [24]. 
This stability equation reduces for n = to the stability 
equation for the emergence of zonal flows [3, 25]. 

For small values of the energy input rate, the growth 
rate 3i(cr) is negative for all structures and the statisti- 
cally homogeneous state is stable and persists. For a crit- 
ical value e c , the homogeneous flow becomes SSST unsta- 
ble, symmetry breaking occurs and coherent structures 



ensemble average represents a Reynolds average with the 
ensemble mean representing coarse-graining. This inter- 
pretation has been adopted in the SSST study of turbu- 
lence in baroclinic flows [22, 23]. With this interpretation 
of the ensemble mean, the SSST system (2)-(3) provides 
the statistical dynamics of the interaction of the ensemble 
average field, which can be a zonal or a non-zonal coher- 
ent structure, with the fine-grained field, represented in 
the theory through its covariance C. The SSST system 
defines an autonomous dynamics and its fixed points de- 
fine a new type of turbulent statistical equilibria. While 
these equilibria formally exist only in the infinite ensem- 
ble limit, it has been shown that their characteristics 
manifest even in single realizations of the turbulent sys- 
tem. The structural stability of these turbulent equilib- 
ria can be addressed in SSST by studying their stability. 
Specifically, when an equilibrium of the SSST equations 
becomes unstable, the turbulent flow bifurcates to a dif- 
ferent attractor. This theory therefore predicts param- 
eters of the physical system which can lead to abrupt 
reorganization of the turbulent flow. 

The SSST equations (2)-(3), admit for v = the equi- 
librium 

U E = V E = 0, C E = =-, (5) 
2r 

that has zero large scale flow and a homogeneous 
eddy field with spatial covariance dictated directly from 
the forcing. We now investigate the SSST stability 
of this equilibrium as a function of the energy input 
rate, e, and relate the outcome of this stability analy- 
sis to the results in the non- linear simulations of (1). 
The stability of the homogeneous equilibrium (5) is as- 
sessed by introducing small perturbations of the form 
[5Z,5C] = [SZ nm ,5C nm ]e m ^ 1+x ^/ 2+ " n ^ 1+ y 2 ^ 2 e' 7t to 
the SSST equations (2)- (3) linearized about the equi- 
librium (5) and calculating the eigenvalue a. When 
3?(ct) > 0, the structure with (x, y) wavenumbers (n, to) 
is unstable and will emerge. It can be shown that the 
eigenvalue a satisfies the equation 
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emerge. The growth rates as a function of the wavenum- 
bers of the structure (n, to) are shown in Fig. 4. For 
e/e c = 2.6, the structure with the largest growth rate, is 
non-zonal with (|n|,|m|) = (1,5) and has 9(cr) = 0.4, 
implying retrograde propagation of the eigenstructure. 
Note also that for this energy input rate, the zonal flows 
are SSST stable and are not expected to emerge. For 
e/e c = 30, both zonal jets and non-zonal structures are 
unstable, but the zonal jets have smaller growth rates 
compared to the non-zonal structures [26]. The zonal 
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FIG. 4. The growth rate, 5R(u) as a function of the wavenum- 
bers \n\, \m\ of the coherent structures for (a) e/e c = 2.6 and 
(b) e/e c = 30 (only positive values of 5R(cr) are shown). Un- 
stable waves with \n\ > have ^s(a) > corresponding to 
retrograde propagation. 




FIG. 5. (a) Snapshot of the streamfunction ip(x,y,i) and 
(b) Hovmoller diagram of ijj(x,y = 7r/4, t) obtained from the 
SSST integrations for e/e c = 2.6. The thick dashed lines in 
(b) show the phase speed obtained from the stability equation 
(6). 

jets have 5s (a) = 0, implying that they grow in place, in 
contrast to the non-zonal coherent structures that always 
propagate in the retrograde direction. Numerical integra- 
tion of the SSST system (2)-(3), shows that for e > e c 
the unstable structures typically equilibrate at finite am- 
plitude after an initial period of exponential growth. As 
a result (2)-(3) admit in general multiple equilibria. Fig- 
ure 5(a) shows the equilibrium structure with the largest 
domain of attraction, when e/e c = 2.6. This structure 
coincides in this case with the finite amplitude equilib- 
rium of the fastest growing (|n|, |m|) = (1,5) eigenfunc- 
tion and propagates as illustrated in Fig. 5(b) in the ret- 
rograde direction with a speed approximately equal to 
the phase speed of this unstable cigenstructurc. A proxy 
for the amplitude of these equilibrated structures are the 
zmf and nzmf indices calculated for the SSST integra- 
tions that are shown in Fig. 1. As the energy input rate 
increases, the non-zonal structures equilibrate at larger 
amplitudes. However, for e > e n ;, the equilibria with the 
largest domain of attraction are zonal jets and the flow 
is dominated by these structures (cf. Fig. 1). 

The results of the SSST analysis are now compared 
to the non-linear simulations. First of all, the SSST sta- 
bility analysis accurately predicts the critical e c for the 
emergence of non-zonal structures in the non-linear simu- 
lations as shown in Fig. 1. The finite amplitude equilibria 
obtained when e > e c also correspond to the dominant 



structures in the non-linear simulations. For e/e c = 2.6, 
the spectra in the nonlinear simulations show significant 
power at (\n\, |m|) = (1,5), the structure predicted by 
SSST to have the largest domain of attraction. Remark- 
ably, the phase speed of the (1,5) waves observed in the 
non-linear simulations and the amplitude of these struc- 
tures as illustrated by the nzmf index are approximately 
equal to the phase speed and amplitude of the correspond- 
ing SSST equilibrium structure (cf. Figs. 1 and 3). For 
e > e n /, in both nonlinear and SSST simulations zonal 
jets emerge and the power of the non-zonal structures is 
substantially reduced. Comparison of the number of jets 
and their amplitude between the SSST and the nonlin- 
ear simulations also shows good agreement. This demon- 
strates that the SSST system can predict the amplitude 
and characteristics of both the non-zonal and the zonal 
structures that emerge in the turbulent flow. 

While the regime transition that occurs at e c is pre- 
dicted by the stability equation (6), the second transition, 
which is associated with the emergence of zonal flows and 
occurs at e n i, is more intriguing. The stability equation 
(6) predicts that the zonal structures become unstable at 
(-sz = 4e c < e„/. In previous studies of SSST dynamics 
restricted to the interaction between zonal flows and tur- 
bulence, these initially unstable structures were found to 
equilibrate at finite amplitude [3, 20]. Preliminary calcu- 
lations show that within the context of this generalized 
SSST analysis that takes into account the dynamics of 
non-zonal structures as well, these equilibria are found to 
be saddles that are stable to zonal but unstable to non- 
zonal perturbations. The threshold for the emergence of 
jets in the SSST simulations and in the nonlinear simu- 
lations is therefore determined as the energy input rate 
at which an SSST stable, finite amplitude zonal jet equi- 
librium exists. A method to correctly obtain the critical 
input rate e n i has been recently developed [20]. It starts 
by recognizing that for e c < e < e sz , the spectrum in the 
turbulent flow is modified (C E ^ S/2r) and is given by 
the covariance C E associated with the finite amplitude 
equilibria similar to the ones shown in Fig. 5. If this mod- 
ification is taken into account, then the stability analysis 
around the equilibrium [Ue,C e ] correctly predicts e„/. 

In summary, SSST constitutes a theory that predicts 
the two regime transitions in the turbulent flow as the 
energy input rate is increased: the symmetry breaking of 
the homogeneous equilibrium with the emergence of co- 
herent, propagating non-zonal structures and the emer- 
gence of zonal jets. It also predicts the characteristics 
of the emerging structures (their scales and their phase 
speed), as well as their amplitude. The relation of the 
non-zonal coherent structures to westward propagating 
vortex rings in the ocean and coherent vortices in plan- 
etary atmospheres will be the subject of future research. 
The fact though that these structures are more suscepti- 
ble to instabilities compared to zonal jets, shows that the 
prevalent structures in planetary flows are zonal jets. 
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